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Abstract 

A phenomenological model for a measurement of "barrier traversal times" for 
particles is proposed. Two idealized detectors for passage and arrival provide 
entrance and exit times for the barrier traversal. The averaged traversal time 
is computed over the ensemble of particles detected twice, before and after 
the barrier. The "Hartman effect" can still be found when passage detectors 
that conserve the momentum distribution of the incident packet are used. 
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The temporal characterization of quantum mechanical tunnelling traces 
back to early studies by McColl 0. More recently a paper by Biittiker 
and Landauer |2| and interest in the subject from various fields (as var- 
ied as nuclear and molecular physics, cosmology or semiconductor physics) 
have triggered a debate that has frequently dealt with the very foundations 
of quantum mechanics ||. The interpretation of the quantum mechanical 
formalism and the wave-particle duality, the quantization ambiguities, the 
relation between classical and quantum mechanics, or the quantum "mea- 
surement problem" are some of the ingredients of this research. These are all 
difficult and not completely understood matters so, not surprisingly, answer- 
ing the question "How long does it take to cross a barrier?", i.e., defining a 
quantum traversal time has been controversial. (We shall mainly discuss the 
general concept of "traversal time" instead of a more restrictive "tunelling 
time". However especial attention will be paid to tunnelling conditions in 
the calculations.) 

In this problem the standard quantization procedures are difficult to ap- 
ply since only a limited number of classical trajectories cross the selected 
region and "continuous observation" may be required for a measurement. 
Even so, many proposals exist that generalize in different formal or opera- 
tional ways the classical concept of traversal time to the quantum case. The 
debate on the barrier traversal time is essentially a consequence of different 
conditionings and criteria, added to the bare original question, that priv- 
iledge one quantum quantity versus the others. As long as the conditioning 
is made explicit, to make clear that different versions of the original ques- 
tion are being answered, there is no fundamental conflict among seemingly 
irreconcilable proposals. (Part of the theoretical work -using path integrals 
[|j] or a projector approach ||- has been devoted to develop comprehensive 
formalisms that allow to classify and relate many of the possible character- 
istic quantities.) However not all aspects have yet been investigated. Only 
the totality of conditionings or additional specificications exhausts the possi- 
ble information about the barrier traversal in the temporal domain. Within 
this spirit we shall investigate here a complementary aspect to those we have 
previously examined [§-10], and to experiments performed on electromag- 



netic waves to measure "Larmor times" fTT|| . The objective of this letter 
is to examine one operational definition of barrier traversal time for parti- 
cles. By "operational" we mean "related to a specific experiment, possibly 
a "Gedanken" experiment. We shall model an idealized experimental setup 
inspired by an elementary "classical" receipe: In order to measure the transit 
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time trough a spatial region the first entrance t\ and first exit times £2 are 
measured and their difference r = t 2 — ti is evaluated for each particle. In 
our case the spatial region includes a potential barrier and only the particles 
detected before and after the barrier will be taken into account. If the exper- 
iment is repeated many times r can be averaged and its statistial properties 
examined. In a previous publication by Muga, Brouard and Sala || a related 
approach was proposed for the quantum case: An average entrance instant 
(t) m at a and an average exit instant (t) out at b were defined in terms of 
incident and outgoing current densities, 

/,\in Jq c J{a)tdt 
/, \ out f °° J(b)tdt 

{t)b - ss°j(b)tdt' [ ) 

with a traversal time tt = (t)° ut — (t)™ given by the difference between the 
two averages. Here J is the current density and it is assumed that a is far 
from the barrier so that the packet passes rightwards through point a before 
t c , a time prior to the backwards reflected flow after the collision. is in 
principle measurable but it has a clear drawback since it is not the average 
of transit times for individual particles. It is instead the difference between 
two averages of different nature. This is better understood in classical terms: 
The average entrance instant (t) m is operationally defined for the ensemble 
of particles that arrive at the first detector while the exit time is only defined 
for a smaller set (those that arrive at the final detector). This definition in 
fact may lead to negative values of in the classical and quantum cases 
[0,0]. Classically the average entrance time may be dominated by trajec- 



tories that are eventually reflected so that (t) m can be very different from 
typical entrance times of the trajectories that eventually pass the barrier. In 
this letter this inconsistency with the classical limit is overcome by restricting 
the averaging to those particles that are detected before and after the barrier. 
In general this approach implies a "back reaction" of the first detector that 
modifies the state. We accept this perturbation as a fact and investigate the 
outcome of the described operational procedure, and the effect of different 
detectors, in particular of those that minimize the back reaction so that the 
momentum distribution of the initial packet is preserved. 

In general the particle+detector system involves many degrees of freedom 
and it is rarely modelled accurately. The objective of a phenomenological 
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model is to retain its essential aspects with the aid of some adjustable set 
of parameters and in agreement with experiental facts. Our model does not 
specify the particular features of the detection at a detailed experimental 
level but we have in mind particle tracks similar to the ones produced in a 
bubble chamber or by means of photographic plates. These tracks are char- 
acterized by a discrete set of macroscopic spots (two in our case) originated 
at certain times ("clicks") considered as "classical events" that result from 
the quantum particle passage or arrival. The particle is restricted to one 
dimensional spatial motion. Specifically the effect of the detector associated 
with a given spot is simulated according to a track formation model proposed 
by A. Jadczyk and Ph. Blanchard using two basic elements |13| : An effective 
one-degree-of-freedom Hamiltonian and a modified projection postulate for 
the particle state after the first detection. 



1 Model description 

The initial state of the particle is given by a wave function ijj associated with a 
preparation procedure. In operational terms, an ensemble of noninteracting 
particles, represented symbolically as {E }, is sent towards the barrier - 
one particle at a time- from the left with identical specifications. (In our 
calculations the initial state at t = is a minimum-uncertainty-product 
Gaussian centered at position x = 20, momentum p = 8 and spatial variance 
9/4, all quantities in atomic units. The potential barrier is a square barrier 
with "height" Vq = 50 from x = 80 to x = 80 + d and the particle has mass 
m = 1.) 

Two particle detectors A and B are located on both sides of a barrier 
potential, at x = a and x = b. The first one is a passage detector that 
does not destroy the particle. The second one is an arrival detector. The 
translational degree of freedom of the particle, x, is the only one represented 
explicitly. A simplifying assumption is that only one of the two detectors 
is working at a time: When the particle is sent to the barrier only A is 
active. Detection of the particle at A disconects this detector and activates 
the second, B. 



4 



1.1 First detector: Probability of detection 

It can be proved using multichannel scattering theory techniques that the in- 
cident channel amplitude (corresponding to translational motion of the par- 
ticle and the detector A in its lower state) can be represented by an effective 
Schrodinger equation with a complex potential fl5fl . (In "Event Enhanced 



Quantum Theory" as described in ]T3| the imaginary part of the potential 



is deduced rigorously from the Lindblad form of the Liouville equation that 
describes a coupling of the quantum system with a classical detector.) Here 
the effective Schrodinger equation is written as 

h 2 d 2 

Hip(x, t) = -——ip(x, t) + [V(x) + A(x)]i/>(x, t), (3) 

where V(x) represents the potential barrier and the complex potential, A, is 
written as 

A{x) = -^g 2 {x;a), (4) 

with 

g{x,a) =se-^~ a)2/2a2 . (5) 

The "intensity" , s, and "width" , a, of the detector are adjustable parameters. 
The norm of the incident channel, 

/oo 
ip*(x,t)if;(x 1 t)dx 1 (6) 
-oo 

decreases, due to the detector presence, from the initial value N(0) = 1. 
The total absorption 1 — iV(oo) is the efficiency of the detector. It is not 
necessarily equal to one so the ensemble of particles detected at A, {E a }, is 
generally smaller than {Eq}. The normalized probability density for trigger- 
ing the detector at time t a is proportional to the absorption rate —dN/dt\t a . 
Normalizing with respect to the ensemble {E a } it is given by 

p , f | cm dN(t a )/dt a 
[ al a) S™dtdN{t)/df [ ) 

1.2 Effect of detection on the wave functions 

It will be assumed, within the spirit of a simplified phenomenological model, 
that after each detection (a "click") the state of the particle can be effectively 
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represented by a modified wavefunction. The true final states should be 
determined by a detailed analysis of the interaction between the system and 
the detector. Instead we shall later assume a physically motivated functional 
form. The ensemble of detected particles can be represented by a statistical 
mixture of such states. This is of course reminiscent of Von Neumman's 
projection postulate. However an important feature of a bubble chamber 
track is that it does not look like a random walk. This cannot be explained 
with a naive projection localizing the particle position by means of position 
eigenstates, since a position eigenstate has equal probability to expand in any 
direction (erasing the memory of the state previous to the measurement) so 
there would be no tendency to ionize atoms in the direction of the dominant 
incident momentum A modified projection postulate correcting this fact 



has been derived by Jadczyk and Blanchard. The wave function resulting 
from a click at time t a and consistent with track formation has a memory 
of the previous state and reflects also the detector properties. A simple 
expression satisfying these two conditions is |L3[ 



, / x g(x)jj(x,t a ) 

where ip(x,t a ) is the wave function evolved with the Schrodinger equation 

To determine the effect of the detector we have examined the momentum 
average and its variance for the ensembles {Eq} and {E a }. Averages over 
{Eg} require some care since they imply a a double average: The first one 
(represented as Q) is a quantum mechanical average using each wave packet 
ip ta ; the second (D) is an average over the times of detection t a weighted by 

P{ta\E a ), 

(p) Ea = DQp = J P(t a \E a )(ij ta \p\ij ta ) dt a . (9) 

Since there are two types of average different "variances" are possible [0,0 • 
For the ensemble {E a } the important one is A 2 D q = DQ[p 2 — (DQp) 2 ]. (This 
is a variance computed over detected particles regardless of their detection 
time [0.) The average momentum is conserved well (especially by weak 
detectors) except for very narrow detector widths. For all detectors used in 
this work DQp « (p)e better than 0.2%. However the "momentum widths" 
A D q (square root of variance) may change drastically with respect to the 
momentum width A p of the original packet. Fig. 1 shows that wider detectors 
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tend to keep the variance of the original state while narrow detectors give 
very large variances. Weak detectors (small s) conserve the variance better 
than strong detectors (large s). In summary, in our model weak and wide 
detectors are the best as far as conservation of the momentum distribution 
of the original packet is concerned. They are however not very efficient, for 
s = 1 the absorbed norm goes from 0.05 to 0.6 in the a-interval of Fig. 1. In 
comparison the full norm is absorbed for s — 10. 



1.3 The second (arrival) detector 

The second detector is assumed to be a perfect one as described in 0, so 
that the full transmitted packet is absorbed. It is located at the right edge 
of the barrier. Let {E b } be the ensemble of particles that produce two clicks 
at times t a and t b and P(E b \t a ) the transmittance of ip ta , i.e., the fraction of 
the norm of tp ta that will be transmitted and therefore detected at B jBJ. 
The probability for being detected at B conditioned to having been detected 
at A is 

P(E b \E a ) = J P(E b \t a )P(t a \E a )dt a . (10) 

Instead of using an expression similar to (^) the distribution of arrival times tb 
for a perfect absorber can be approximated accurately by the (normalized) 
flux without absorber |J. In particular, for a wave packet ifj ta {x',t a ), the 
detection probability density at t b in B, conditioned to having been detected 
at t a in A and restricted to the ensemble {E b }, is given by 

WtJ = TOTS' (u) 

where J ta is the flux for the state ipt a - Using Bayes' rule the joint probability 
density for detection at t a in A and t b at B restricted to the ensemble {E b } 
is given by 

p( . , ipx P(t b \E b ,t a )P(E b \t a )P(t a \E a ) 

P(t6 ' talEb) ~ I P(E b \t a )P(t a \E a )dt a ■ (12) 

Finally, the probability distribution of r = t b — t a is computed, for the 
ensemble {E b }, by integrating over t b and t a with the delta function 5{t b — 
ta - t), 

p( r |£ 6 ) = [ P ( ta ± T \ Eb ' tg)P(E b \t a )P(t a \E a )dt a 
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We have calculated average traversal times (r) E = / P{r\E^)rdT versus the 
barrier width d for two different weak detectors at a, both with s = 1. One 
of them, Ai, is a wide one and conserves well the momentum distribution of 
{E }. The other one, A 2 , is a narrow detector, and produces a momentum 
variance which is approximately ten times the initial one. The detector before 
the barrier is always put far from the barrier (a = 50) to compare with 
the type of Gedanken experiment performed in ref. ||, so that the initial 
packet may pass through a before interacting significantly with the barrier, 
and b is located at the right barrier edge. Let T\ and r 2 be the averages 
corresponding to using the two initial detectors A\ and A 2 . Figure 2 shows 
that the Hartman effect, i.e. the fact that the average traversal time does not 
grow with d (actually it decreases slowly [0]) can still be seen with A x until 
a critical barrier width d c where the "classical passage" of momenta "above" 
the barrier starts to dominate [||,|| . When the narrow detector A 2 is used the 
momentum variance is so large that the transmission is always dominated by 
fast momenta well above the barrier (We have independently checked this 
fact by calculating the ratio between transmission due to energies above 
and below the barrier energy.) so that the behaviour is the one expected 
classically, i.e., a linear growth of t 2 with d. Figure 2 also shows tt, which 
is qualitatively very similar to t\. The relation tt < T\ is due to the two 
different ways the average is performed in the initial detector and can be 
also understood on classical grounds. The right front of the incident packet 
is dominated by faster momenta and it contributes with more particles to 
the transmitted ensemble. For computing the later, no distinction is made at 
a between particles to be transmitted or not. (The effect grows with d until 
it saturates when the transmission is purely above the barrier.) Note that 
tt could be negative while the times defined in the present work are always, 
by construction, strictly positive. The "displacement" of the curve r a with 
respect to tt (note the difference in the value of the critical barrier) is due 
to the slight difference in the momentum variances. 

In summary a two detector measurement of a particle traversal time has 
been modelled. Passage detectors conserving the initial wave packet momen- 
tum distribution still show the Hartman effect. 
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Figure Captions 

Figure 1. Square root of the momentum variance after detection, A D q, 
for s = 1 (solid line) and s = 10 (dashed line). The dashed-dotted line is the 
reference value of the momentum variance for the original ensemble {E }. 

Figure 2. Average traversal times versus barrier width d evaluated for 
(a) s = 1, a = 4.5 (dashed line); (b) s = 1, a = 0.2 (dashed-dotted line). 
The average time tt is also represented (solid line). 



11 



